ANALYSIS OF A FORCE-BASED QUASICONTINUUM APPROXIMATION 



MATTHEW DOBSON AND MITCHELL LUSKIN 

Abstract. We analyze a force-based quasicontinuum approximation to a one-dimensional system 
of atoms that interact by a classical atomistic potential. This force-based quasicontinuum approx- 
imation can be derived as the modification of an energy-based quasicontinuum approximation by 
the addition of nonconservative forces to correct nonphysical "ghost" forces that occur in the atom- 
istic to continuum interface. The algorithmic simplicity and improved accuracy of the force-based 
quasicontinuum approximation has made it popular for large-scale quasicontinuum computations. 

We prove that the force-based quasicontinuum equations have a unique solution when the mag- 
nitude of the external forces satisfy explicit bounds. For Lennard- Jones next-nearest-neighbor 
interactions, we show that unique solutions exist for external forces that extend the system nearly 
to its tensile limit. 

We give an analysis of the convergence of the ghost force iteration method to solve the equilibrium 
equations for the force-based quasicontinuum approximation. We show that the ghost force iteration 
is a contraction and give an analysis for its convergence rate. 



1. Introduction 

The local lattice structure for minimum energy configurations of atomistic systems subject to 
external forces is usually slowly varying except near defects such as dislocations [12]. Quasicon- 
tinuum methods efficiently approximate these multiscale features by maintaining atomistic degrees 
of freedom near defects and coarse-graining the atomistic degrees of freedom in regions where the 
local lattice structure is nearly uniform through the introduction of representative atoms [7,12,23]. 
The efficiency of quasicontinuum methods has allowed the simulation of more complex problems 
than can be computed using a completely atomistic model [21]. 

Many quasicontinuum methods have been proposed [5,7,11,12,17,19,22,23], and each version 
gives a different quasicontinuum approximation of the atomistic system. A force-based quasicontin- 
uum approximation has been proposed that modifies an energy-based quasicontinuum approxima- 
tion by the addition of nonconservative forces to correct nonphysical "ghost" forces that occur in 
the atomistic to continuum interface [7,12,19,21]. The force-based quasicontinuum approximation 
has been popular for large-scale computations since the improved accuracy is obtained with no 
additional computational work simply by computing the force on each representative atom with 
either an atomistic algorithm or with a continuum finite element algorithm [7,12,18,19,21]. 

Adaptive mesh and error control have been successfully used with the force-based quasicontinuum 
approximation to efficiently choose representative atoms [7,15,18,21]. The number of representative 
atoms surrounding defects that need to be modeled atomistically (the core of the defect) can be 
determined by the error tolerance, and the mesh in the continuum region surrounding the core can 
be coarsened beyond the atomistic-continuum interface wherever the deformation gradient varies 
slowly. For simplicity, reported implementations have not coarsened within the cut-off radius of 
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atomistic representative atoms, but it is possible to coarsen immediately beyond the atomistic- 
continuum interface by interpolating between continuum representative atoms. 

In Section O we give a derivation following [7, 12, 19, 21] of several quasicontinuum approxima- 
tions leading to the derivation of the force-based quasicontinuum approximation. In Section [3l we 
reformulate the equilibrium equations as a balance of forces conjugate to the distances between 
representative atoms, rather than as a balance of forces conjugate to the positions of the represen- 
tative atoms. Our derivation and reformulation gives the mathematical structure that is used in 
our analysis. 

In Section U we prove that the force-based quasicontinuum equations have a unique solution 
under suitable restrictions on the loads. In the case of Lennard-Jones next-nearest-neighbor inter- 
actions, we determine bounds for the magnitude of the loads for which unique solutions exist and 
find that the allowable loads extend quite close to the tensile limit. 

In Section \5\ we give an analysis of the convergence of the ghost force iteration method that has 
been most commonly used to solve the equilibrium equations for the force-based quasicontinuum 
approximation [7,12,19,21]. We prove that the ghost force iteration is a contraction and give a 
bound for its convergence rate. We show that our bound for the convergence rate gives a high 
convergence rate when applied to the Lennard-Jones model subject to moderate external forces. 

Mathematical analyses of energy-based versions of the quasicontinuum approximation that do 
not include ghost force corrections have been given in [2,8,9,13,14,16,17], and a simplified version 
of our analysis can be used to prove the existence of solutions to these energy-based quasicontin- 
uum approximations. We show, though, that the ghost forces are nonconservative forces, so they 
cannot be derived from an energy. Thus, the force-based quasicontinuum approximation cannot be 
completely analyzed by energy methods. 

We refer to [3] for a review of current progress on the mathematical analysis of atomistic to 
continuum models for solids and to [12] for an introduction and overview of the quasicontinuum 
approximation . 



2. Quasicontinuum Approximations 

In this section, we describe a sequence of one-dimensional coarse-grained approximations of a 
chain of atoms with nearest-neighbor and next-nearest-neighbor interactions given by a classical 
two-body potential, (j){r). We assume that the atomistic potential (j){r) is defined for all r > 0. 

We begin with the atomistic model, which has degrees of freedom for all atomic positions and 
computes a total internal energy directly from pairwise interactions. From there, we examine 
the constrained atomistic approximation and the local quasicontinuum approximation which both 
decrease the degrees of freedom by interpolating atomic positions between representative atoms. 
We then introduce an energy-based and a force-based quasicontinuum approximation which span 
atomistic and continuum scales by combining atomistic regions where the atoms directly interact 
according to the atomistic model and continuum regions where the atoms interact according to the 
local quasicontinuum approximation. We observe that the energy-based quasicontinuum approxi- 
mation gives nonphysical ghost forces near the atomistic to continuum interface that are corrected 
by the force-based quasicontinuum approximation. 



y-M Vj-i Vj Vj+i vm+i 



Figure 1 . Atomistic chain with atoms labeled by their position. 
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2.1. The Atomistic Model. We denote the positions of the atoms by yi for i = —M, . . . , M + 1, 

where Ui < yj+i. The total energy for the atomistic system with nearest-neighbor and next-nearest- 
neighbor interactions given by the classical two-body potential (/)(r) is 



M 



s^iy) = [<P{yi+i - Vi) + (t}{yi+2 - vi)] 



i=-M 

where y = {u-m, ■ ■ ■ ,yM+i) G m2A/-i-2 where the boundary terms 4>{yi — yj) above and in the 
following should be understood to be zero for i ^ {— M, . . . , M + 1} or j ^ {—M, . . . ,M + 1}. We 
can also express the total energy in terms of energies associated with each atom as 

M+l 

^'^(y) = E ^^(y) (2.1) 

i=-M 

with 

^iiy) = 2 ^(yi+^ - Vi) + (f'iyi+2 - Vi) + 4>{yi - yi-i) + 4>{yi - ^^-2) • (2.2) 

We then have that the force on the atom at position yi is given by 



d r 

-^"(y) = [Hyi+i - yi) + 4>{yi+2 - yi) + - yi~i) + Hyi - yi-2) 

= [r]{ri) + r]{ri + r^+i)] - [r/(ri_i) + r/(ri„i + ri-2)] , 



(2.3) 



where T]{r) = (j)'{r) and rj = —yi is the lattice spacing at yi. The terms ?7(rj) and r]{ri+rj) above 
and in the following should be understood to be zero for i ^ {—M, . . . , M} or j ^ {— M, . . . , M}. 

We now assume that the atoms are also subject to an external force, fi{yi), that is obtained from 
an external potential energy of the form 

A/+1 

r^iy) = y. vfiy,), 



i=-M 



SO 



For example, such external forces may model the interaction of the one-dimensional chain with 
atoms in layers above and below the chain, as in the Frenkel Kontorova model [10]. 
We then have the equilibrium equations 

i^"(y) + /i(2/i) = 0, i = -M,...,M + l. 

2.2. The Constrained Atomistic Quasicontinuum Approximation. One can reduce the de- 
grees of freedom in the atomistic model by linearly interpolating the positions of the atoms between 
a set of representative atoms (see [17] for an analysis in this case). We introduce representative 
atoms with positions Zj such that 

zj=y,^ fori = -iV,...,A^+l, (2.5) 

where i-N = — M, In+i = M + 1, and < ij+i- We then let Vj = £j+i — £j denote the number of 
atoms between Zj and Zj^i (see Figure [2]). We now have that 
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is the distance separating i/j equally spaced atoms between zj and zj-^-i, and we have the conservation 
of mass equation 

N N 

E = E (^i+i - = ^^+1 - = 2M + 1. (2.6) 

j=-N j=-N 



y-M yij y^j+i VM+i 



Figure 2. A coarsening of the atomistic chain. 



In the constrained atomistic model, the positions of atoms between Zj and Zj+i are linearly 
interpolated as 

u ■ — ^ 2 

yij+i = yij+i-{v,-i) = — — Zj^ — zj+i forO<i<z/j, (2.7) 

and we define the total internal energy in terms of z = (z-at, . . . , ^at+i to again be the 

interaction enegy of all atoms in the chain, computed according to (|2.ip . giving 

M+l 

f^(z) = £:;(y(z)). (2.8) 

j=-M 

So far, we have reduced the degrees of freedom necessary for denoting the atomistic positions, 
but the total energy is still computed as a sum of energy contributions from all of the atomistic 
degrees of freedom. However, all nearest-neighbor and next-nearest-neighbor contributions for 
atoms interpolated between a pair of represenative atoms are identical due to the uniform spacing, 
so if Ui does not denote a representative atom we have that 

^{Vi+i - Vi-i) = (PC^iyi+i - Vi)) = 4>{'^{yi - Vi-i)) = 0(2r-j) for <i < ij+i. 

These interactions account for all of the energy contributions except for next-nearest-neighbor 
interactions that straddle a representative atom. We can treat these interactions by observing that 

1 1 

4>{yi+i - Vi-i) - ^H'^ivi - Vi-i)) - 2<^(2(yi+i - Vi)) 

_ \(p{rj-i +rj) - i</>(2rj_i) - ^(pi2rj), if yi = Zj for some j = -N,.. . ,N+1, 
[0, if / Zj for all j = -N, ...,N + 1. 

We can therefore partition the total energy into the energy of the region between Zj and Zj^i 
for each j = —N, . . . , N, plus interfacial energy terms that account for interactions that straddle a 
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representative atom. We have from ()2.8p that 



M 



M 



"^""(z) = ^ <P{yi+i -yi)+ X] 4>{yi+i - Vi-i) 

i=-M i=-M+l 
M 

= J2 ['^(^'■+1 ~ + <^(2(yi+i - Vi))] 

i=-M 

M 



i=-M+l 



Hvi+i - Vi-i) - ^(PC^ivi - Vi'i)) - ^0(2(yj+i - Vi)) 



(2.9) 



N 



^0(2(y-A/+i - V-m)) - ^(/)(2(yM+i - Vm)) 

N+l 



-N 



-N 



where 
and where 



4>{r) + 4>{2r) for r > 



(2.10) 



5^Af(r_Ar) 



1 

2 

^0(2rjv). 



(2r, 



-iV + l,...,Af, 



We note that (p{r) = <j){r) + 0(2r) is the energy per atom for an infinite atomistic chain with the 
uniform lattice spacing j/j+i — = r for ah -co < i < oo, and that Sj{rj^i,rj) can be considered to 
be a surface energy at Z-n and z^^i and to be an interfacial energy at Zj for j = — 1, . . . ,N. We 
observe that Sj{rj^i,rj) is a second divided difference for 0(r) about r = rj_i + rj with increment 



rj-i, so 
5j(rj_i,rj) 



1 



)"(r,„i + rj){rj - rj.^f + 0{\rj - r^-il^), j 



-7V + i,...,Ar. 



We will see that if (j){r) satisfies the assumptions given in Section HI then the convexity condition 

Sj{rj_i,rj)>0, j = -N+l,...,N, 

holds for the range of r = (r_Ar, . . . , tat) G M^^+i defined by (|i36|) in TheoremOand CorollaryO 
where solutions of the force-based quasicontinuum equilibrium equations are shown to exist and 
where the iterates of the ghost force iteration reside. 

2.3. The Local Quasicontinuum Approximation. If we neglect the surface and interfacial 
energy terms, Sj, in (j2.9p . then we obtain the local quasicontinuum approximation [12] 

Af+l 



j=-N 



where 



(2.11) 
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To treat the boundary terms consistently, we set U-n-i = ^N+i = 1 and </>(r_Ar_i) = (j){r]\f^i) = 0. 

We remark that the representative atoms need not be placed at atomistic sites as described by 
()2.5p . The local quasicontinuum approximation only requires that Zj < Zj+i for j = —N,. . . ,N 



and that the uj are positive and satisfy the conservation of mass condition (|2.6p . The approxima- 
tion can be generalized to higher space dimensions by using the Cauchy-Born rule [12]. The local 
quasicontinuum approximation is computationally simpler than the constrained atomistic quasi- 
continuum approximation, especially in higher dimensions where the computation of the interfacial 
energy becomes expensive. In the following, we will sometimes refer to the local quasicontinuum 
approximation as the continuum approximation. 

The force on a representative atom at zj for j = —N -|- 1, . . . , is given for the local quasicon- 
tinuum approximation by 



d 




dzj 


H 


d 




dzj 


H 




+ i 



Vj + 



Zj-1 



2{zj+i - Zj) 



(2.12) 



Vj-l + 



2{z, 



[r/(r,) + 2ry(2r,)] - [r?(r,-_i) + 277(2r,-_i)] 



where again 



■i+i 



is the lattice constant for the atoms between Zj and 2:^+1. In the above, we see the "local" nature 
of the approximation, as the force on a degree of freedom is determined only by the positions of 
adjacent degrees of freedom and no long-rage interactions occur. 

We can similarly compute the force on the boundary atoms, noting the one-sided nature of 'S'^(z) 
for the boundary atoms, by 



d 



dz^N 
d_ 

dz_N 
[ri{r-N) 



Z-N+l — Z-N 
Z-N+l - Z-N 



F, 



+ 2ri{2r 
d 



dzN+1 

d 

dzN+i 



ZN+l — Zn 
/ ^iV-i-l - Zj^ 



+ 



2(z_ 



N+l 



z~n) 



(2.13) 



UN 



+ 



2{zN+i - Zn) 



UN 



UN 



(2.14) 



= - [v{rN) + 2v{2rN)] • 

We note that the local quasicontinuum energy, £^{z), and the forces, F^^(z), depend only on 



r = {r-N, ■ ■ ■ jTn), and we will denote the dependance by f (r) and Fj(r) without introducing 
distinct functions. 



We can also derive a local quasicontinuum approximation for the external potential, V (z), by 
setting 

V^iyiz)). (2.15) 



r'-iz) 
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By (j2.4p and (|2.7p . the external force on the representative atom at position Zj is 

"3-1 / ■\ ^3 



=0 

It follows from the linear interpolation (|2.7p that 

We shall assume in our analysis that the external forces, /j, are independent of y. In this case, 
the local quasicontinuum forces, /j, are independent of z and 

/. = E f ^) 4-. + E f ^) , = -iV, . . . , iV + 1. (2.16) 



1=0 



We consider any term fj to be zero if j ^ {— , . . . , M + 1} and any term i^j to be one if j ^ 
{-iV,...,iV}. 

2.4. The Energy-Based Quasicontinuum Approximation. To describe the energy-based qua- 
sicontinuum approximation [12], we again introduce representative atoms with positions Zj for 
j = —N, . . . , N + 1, where Zj < zj^i. Each representative atom is considered to be an "atomistic" 
or "continuum" degree of freedom and will contribute either £j{z) or ^j'(z) to the total internal 
energy according to the atomistic model p.2p or the local quasicontinuum approximation ()2.1ip . 
respectively. 

In applications, atomistic degrees of freedom are used in regions of interest where highly non- 
uniform behavior is expected. Continuum regions surround this, gradually coarsening by increasing 
I'j in regions with slowly varying strain. For simplicity of exposition, we will consider an approxima- 
tion with a single atomistic region, symmetrically surrounded by continuum regions large enough 
so that no atomistic degrees of freedom interact with the surface atoms through nearest-neighbor 
or next-nearest-neighbor interactions. 

We denote the representative atom positions by Zj and define the range j = —K + 1, . . . , X to 
be atomistic sites and the ranges j = —N, . . . , —K and K + 1, . . . , N + 1 to be continuum sites. 
Therefore, the total quasicontinuum energy, £^'~''{r), for the chain is given by 

-K K N+l 

^«^(r) = E ^/(r) + E W + E ^/(r)' (2.17) 

j=-N j=-K+l j=K+l 

where £j{Tc) is defined in (j2.2p and £j'{r) is defined in (j2.1ip . We assume that Uj = 1 for j = 
—K — 1, . . . ,K + 1. This guarantees that Vj = 1 within the next-nearest-neighbor cutoff radius of 
any atomistic site and enables a seamless transition to the continuum approximation. 



— o 

Figure 3. One end of the quasicontinuum chain, highlighting the interface. Filled 
circles are atomistic representative atoms, whereas the unfilled circles are continuum 
representative atoms. 
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The force, F^^^{t), at Zj is then given by 



-,QCE 



■' l=-N 



dr£ dz. 



= < 



F^K+2i^) + H^-K+r-K+i), 



F^{r), 



F^(r) + r]{2rK) - \r]{rK + vk+i), 
Fk+iir) + vC^rx) - ^virx-i + tr) 
Fk+2i^) - hi'^K+TK+i), 
FHr), 



-N < j < -K - 2, 
j = -K- 1, 



j = -K, 

= -K + l, 
j = -K + 2, 
-K + 3 < j < K -2, 
j = K-l, 
j = K, 
j = K + l, 
j = K + 2, 
K + 3<j<N + l. 



(2.18) 



In the above expression, we notice that in the large ranges interior to the atomistic and continuum 
regions the forces are exactly those from the individual models, namely, either F^{r) or Fj'^r). Near 
the atomistic to continuum interface, there are additions to these force terms which contain non- 
physical "ghost" forces. 



• • o o 

ZK-1 Zk ZK+1 Zk+2 

Figure 4. Direction of nonzero forces for a uniform configuration. 



To see this, we consider the forces on the representative atoms when the lattice spacings are 
uniform, that is, when rj = a for j = —N, . . . ,N, ot equivalently when r = a = (a, . . . , a) € M^^^^. 
We observe that the forces computed according to the atomistic model (12. 3p or the continuum 
model (j2.12p are zero except at the ends of the chain, that is, F^{aL) = for j = —N + 1, . . . ,N 
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and F°(a) 



for j = -N + 2, . . . , iV - 1. We then have from ([2J8]) that 



f: 



QCE 



(a) 



[7/(a) + 277(2a)] , 
0, 

ir?(2a), 

-ir?(2a), 

-ir/(2a), 

0, 

-ir/(2a), 

i7?(2a), 

ir?(2a), 

-ir/(2a), 

0, 

- [7?(a) + 2r7(2a)] 



3 



-N, 

-N + 1 < j < -K - 2, 

j = -K- 1, 

3 = -K, 

j = -K + 1, 

j = -K + 2, 

-if + 3 < j < - 2, 

3=K-1, 

3=K, 

j = K + l, 

j = K + 2, 

K + 3<j <N, 

j = N + l. 



However, a continuum chain with uniform lattice spacings has forces only at the surfaces. Figured] 
shows the ghost forces at one of the interfaces. Figure [5] shows the origin of out of balance forces for 
a single representative atom. We note that while the first neighbor terms are balanced, the second 
neighbor terms are not balanced due to the fact that the continuum site does not contribute any 
second-neighbor interactions. 



r?{2a) ri{a) 



ir,{2a) 



Figure 5. Imbalance of forces on an atomistic representative atom near the interface. 



The force-based quasicontinuum approximation described in the next subsection corrects these 
ghost forces. 

2.5. The Force-Based Quasicontinuum Approximation. The force-based quasicontinuum 
approximation [12] corrects the nonphysical forces described in the previous subsection. The forces 
on the representative atoms in the interior of the atomistic and continuum regions are defined by 
(I2.3P and (I2.12P , but the forces on the representative atoms near the atomistic-continuum transition 
given by the energy-based quasicontinuum method must be modified to remove the non-physical 
terms. 

In the force-based quasicontinuum approximation, we again partition into atomistic and con- 
tinuum representative atoms, where the force on a representative atom is the force that would 
result on it if the approximation was entirely of its respective type (atomistic or continuum). With 
this convention, a continuum representative atom only interacts with adjacent degrees of freedom 
regardless of how close any atomistic sites may be. We will see that the tradeoff for this simple 
philosophy is that the forces are not conservative, that is, they cannot be derived from an energy. 

We now model the forces on the representative atoms for j = —K + 1,...,K where K < N — 1 
by the atomistic model (j2.3p and the forces on the representative atoms for j = —N, . . . , —K and 
j = K + l,...,N + lhy the local quasicontinuum approximation (j2.12p - (j2.14p . The forces on the 
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representative atoms are then given by 



F, 



QCF 



(r) 



if(r), -N<j<-K, 
F/(r), -K + l<j<K, 
F/(r), K + l<j<iV + l, 

[v{r-N) + 2??(2r_jv)] , 

[r/(r,) + 277(2r,)] - [virj-i) + 2r/(2r,_i)] , 

[r]{rj) + i]{rj + r^+i)] - [nirj-i) + virj-i + rj-2)] , 

[riirj) + 2n{2rj)] - [r?(r,_i) + 2r?(2r,_i)] , 

- [ri{rN) + 2r7(2rjv)] , 



J = -N, 

-N + l<j< -K, 
-K + l<j<K, 
K + l<j<N, 
j = N + l. 



The force-based quasicontinuum formulation has the desired property that if we take a uniform 
configuration, a = (a, . . . , a) G ]R^-^+^, we have that 



Ff^^(a) = 0, j = -N+l,...,N, 



and on the boundary we get values of equal magnitude and opposite signs. Thus, the equilibrium 
equations have a uniform solution, a, whenever the external forces are applied at the boundary, 
that is, a chain in uniform tension or compression. 

However, the solution to the equilibrium equations for the force-based quasicontinuum method, 



F^^^{r) + fj = 0, j = -N,...,N + l, 



(2.19) 



cannot be obtained from the minimization of an energy since F^^^ {y) is a nonconservative force. 
To see this, we observe that the forces given by the force-based quasicontinuum approximation, 
(r), are not the differential of an energy, since F- (r) is not a closed form [6]. 

We can see this by noting that 



dzK+i 
dzK 



(r) =r?'(rif)+4r/'(2rif). 



implies that 



fiwQCF QpQCF 



dzK 



2.6. The Ghost Force Iteration. In the quasicontinuum method of [12], the forces. F^'^'^{r), 
are split into the force from the energy-based quasicontinuum approximation, F'^'-^^(r), and ghost 
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force corrections defined by 



■f: 



QCE 



(r) 



0, 

+r]{2r^K) - ^V{r-K + r-K+i) 
+?7(2r_x) - ^r]{r^K-i + r-x) 

0, 

+ ir/(rx-i + rx), 
-r?(2rx) + ^-nirK + rx+i), 
-??(2rx) + \-r]{rK-i + rx), 
+ ^'n{rK + rx+i), 
0, 



-N < j < -K - 2, 
j = -K- 1, 
J = -K, 
j = -K + 1, 
j = -K + 2, 

-K + 3 < j < K -2, (2.20) 

j = K-l, 

j = K, 

j = K + l, 

j = K + 2, 

K + 3 < j < N + 1. 



The forces Fj^{r) act as a model correction near the atomistic-continuum interfaces to enforce the 
convention that each representative atom has forces acting on it as though it were surrounded by 
representative atoms of the same type. 

In [12], the solution to the equilibrium equations ()2.19p is obtained by solving the iteration 



(r"+i) + i^^(r") + /,=0, 



(2.21) 



by using a conjugate gradient method to compute r"+^ from r". If the sequence of solutions {r"} 
converges, then the iterative limit r satisfies the equilibrium equations (I2.19p . 

We have split the internal forces of the force-based quasicontinuum approximation, F^'^'^{r), into 
a conservative force, F^^^ {y), and correction, F^{y), to define the iterative method above ()2.2ip . 
The conservative force, F^^^{t), was defined by the quasicontinuum energy, f'^'^(r), given by 
()2.17p . and the correction, Fj^(r), was defined simply by Fj^^(r) = Fj^'^^(r) — F^^^'^{r). However, 
we are free to try to improve the rate of convergence of r" ^ r by constructing an energy £^'-'^'^[r) 
different from £^^^{r) and then solving the iteration 



'(r"+^) + ^(r") + /, = 0, j = -N,...,N + l, 
where F^^^'^{r) is now the force that is derived from £^'^^'^[r) and Fj^^(r) = f'^^^ { 



-f: 



QCE2 



(r). 



3. Conjugate Forces 

In this section, we simplify the analysis by reformulating the force-based quasicontinuum equi- 
librium equations (j2.19p in terms of forces conjugate to the distance between representative atoms, 

Rj = Zj+i - Zj = UjVj, j = -N, ...,N, 

rather than conjugate to the positions of the representative atoms, zj. We will use the notation 
R = {R-N, ■ ■ ■,Rn) G M2^+^. Through this technique, we will be able to derive equations that 
are decoupled inside the continuum regions and are the sum of tridiagonal terms and nonlocal 
interfacial terms in the atomistic region. 
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3.1. The Internal Conjugate Force. We define the internal conjugate force for the energy-based 
quasicontinuum approximation by 



(r) = (r(R)) = -^r-(r), 



dRi 



-N,...,N. 



(3.1) 



We note that we have found it convenient to define i/^f{v) as the negative of the usual convention 
for a conjugate force. We next derive the following relation between F^^^{r) and the internal 
conjugate force ipf{r) : 

QgQC 



f; 



QCE 



(r) 



dzj 



i'fir] 



where we set 



-(r) 

V'f-i(r) 



^r(R)) ' 



dR, 



dzj dRj^i 
N,...,N + 1, 



(r(R)) ' 



dzj 



(3.2) 



L(r)=^f+i(r)=0. 



_7v-ilrj 

We can sum the forces from the left of the chain and use the preceding equation 
that 

(r) = E 

i=-N 

We can derive from either (13.11) or (13.31) that 



J 



-iV, 



,N + 1. 



to obtain 



(3.3) 



<(r) 



'o, 




j = -AT - 1 


Virj) 


+ 2r?(2rj), 


-A^ < J < -i^ - 


Virj) 




j = -K- 1, 




+ Ivirj + rj-i) + \r]{rj + r^+i) + r/(2rj), 




V{rj) 


+ ^r/(rj + rj_i) + r/(rj + r^'+i), 


j = -K + 1, 


< ri{rj) 


+ r?(rj + rj_i) + ??(rj + r^+i), 


-iT + 2 < j < 


virj) 


+ 77(rj + rj_i) + ir/(rj + r^+i), 


j = - 1, 


virj) 


+ ^r/(rj + rj_i) + \r]{rj + r^+i) + r/(2rj), 




virj) 


+ 2r/(2rj) + |r/(rj+rj_i), 


j = K + l, 


virj) 


+ 2r?(2r,), 


K + 2 < j < N, 






j = N + l. 



(3.4) 



We cannot properly derive an internal conjugate force for the force-based quasicontinuum ap- 
proximation since it is not a conservative force. However, we will find it convenient to define 
an internal conjugate force for the force-based quasicontinuum approximation by following p.3p 
setting 





and 



^^N-iir) 



QCF, 



j = -N,...,N +1. 



(3.5) 



-N 



We can then obtain that 



7QCF, 



(r)=Vf(r)-<_i(r), 



-iV,...,iV + l. 



(3.6) 
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We have the fohowing closed form expressions for the internal conjugate force ^J'(r) : 



<(r) 



fo, 

V{rj) + 2r/(2rj), 

r]{rj) + r]{rj + rj^i) + r]{rj + rj+i) 

+ [2r?(2r_x) - vir-K + T-k-i) - vir^K + r^x+i)], 
r]{rj) + 2r]{2rj) 

+ [2r/(2r_ii-) - r]{r_K + T-k-i) - r]{r^K + r-K+i)] 

-[2r?(2rx) - + rx-i) - r]{rK + r-^+i)], 
[2r]{2r^K) - r]{r^K + T-k-i) - r]{r^K + r-K+i)] 

-[2r]{2rK) - r]{rK + rx-i) - r]{rK + tk+i)], 



j = -N- 1, 
-N <j < -K, 

-K + 1 < j < K, 



K + 1 < j < N, 



j = N + l. 



(3.7) 



The internal conjugate force, ^^(r), takes a simpler form when there is no resultant quasicon- 



tinuum force, that is, when r satisfies 



AT+l 



(3i 



j=-N 

which is equivalent to 

2r]{2r_K) - r]{r-K + r^K-i) - r]{r^K + r-K+i) 

= 2r]{2rK) - r]{rK + rx-i) - r]{rK + rx+i). 

We note that (13. 9|) is satisfied when r is symmetric, that is, when r_j = rj for j = 1, . . . ,N. Let 
us now define V^J' (r) on all of IR^^+^ as a symmetric extension of V'j^(r), with equality whenever r 
satisfies (13. 9p . This leads to ipf{r) having a more symmetric form. 



{3S 



0, 



j = -N- 1, 
-iV < j < -K, 



(3.10) 



??(rj) + 2r]{2rj), 

7]{rj) + 7/(rj + rj_i) + 7]{rj + r^+i) 

+ [2r/(2r/^) - ri{rK + rK-i) - ri{rK + tk+i)], -K+l<j<K-l, 
7?(r,)+2r/(2rj), K < j < N, 

0, j = N + l. 

The intervals of definition changed slightly from (j3.7p , which is one of the simplifications afforded 
by (j3.9p . We use ipf{r) in our subsequent analysis, and in Section [331 we discuss the relation 
between using ^pj'{r) and ijjj {v) to solve the equilibrium equations (j2.19p . 

We identify in (|3.10p a continuum internal conjugate force for j 
given by 

r/(rj) + 2ri{2rj) 

and an atomistic internal conjugate force for j = —K + l,...,K — l given by 

r]{rj) + r/(rj_i + rj) + r]{rj + rj+i). 
We identify the remaining terms, 

2r]{2rK) - r]{rK + tk-i) - r]{rK + rx+i), 
for j = —K + 1, . . . , X — 1 as the nonlocal part of the internal conjugate force. 



-N,...,-K and j = K,...,N 
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For consistency, we wish to define ipj{r) as we did for ^j^(r). Since V'j (r) is derived from the 
energy £^'3'^^(r), it has no resultant force. Thus, we define = ipj'{r). 

We can also derive a corresponding internal conjugate force, V'j'(r), by summing Fj^(r) as in (j3.5p . 
From the definition of Fj^(r) ([2:20]) . we have that 



^f(r) = <(r)+Vf(r), 



-N 



We can define ij^fir) by 

V;f(r) = ^f(r)+Vif(r), j 
and we can then check that 
'0, 

-\r]{r_K + r-K-i), 

ri{2r^K) - ^Vir-K + r-K+i) - Ivir-K + T-k-i 

2T]{2rK) - v{fK + TK-l) - VifK + TK+l), 

2r]{2rK) - \'q{rK + rx-i) - iJ^rK + tk+i), 
r]{2rK) - \'n{rK + rx-i) - ^rji^K + rx+i), 
-^vi^K +rK+i), 
0, 



iV-l,...,iV + l. 
1,...,A^ + 1, 



-N -1 < j < -K -'z 

j = -K- 1, 

3 = -K, 

j = -K + 1, 

-K + 2 < j < K -2, 

j = K-l, 

j = K, 

j = K + l, 

K + 2 < j < N + 1. 



(3.11) 



If r satisfies the condition (j3.9p . then we have that ipf{r) = ^^(r). We note that the nonlocal part 
of the internal conjugate force is found in ipj'ir) and ■ip^{r). 

We finally observe that if r is symmetric, then V'j^(r) = ■0j^(r), ■(/'^'(r) = ■ipj'lr), and '(/'^(r) = 
ip'j{r) are symmetric. 

3.2. The External Conjugate Force. We recall that we are assuming in our analysis that the 
external forces, /j, are independent of y, and that consequently the local quasicontinuum forces, 
/j, given in (I2.16P are independent of z. The external potential, ■p^(z), given in ()2.15p thus has the 
form 

Af+l 
j=~N 

We now assume that there is also no resultant force from the external forces, so that 

TV+l 

E = 0- (3.12) 

j=-N 

It then follows that the external potential, V^{z), is also a function of r, and we can define the 
external conjugate force by 



(r(R)) 



1 dV^ 



(r), 



-N,...,N. 



We next derive the following relation between fj and the external conjugate force $j : 



dzj 



dR, 



(r(R)) 



(3.13) 



(3.14) 



J 



-TV, ...,A^ + 1, 
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where we set 

$_7V_1 = ^>7v+i = 0. 

We can sum the external forces from the left of the chain and use the preceding equation (|3.14p to 
obtain that 

i 

= - j = -N,...,N+ 1. (3.15) 

i=-N 

If the external forces, fj, are anti-symmetric, that is, 

/,+i = -/_„ j = 0,...,N, (3.16) 

then we can conclude from (|3.15p that the external conjugate force, <I>j, is symmetric about j = 0, 
that is, 

= j = -N-l,...,N + l. (3.17) 

We note that the external forces, fj, are anti-symmetric if the chain is subject only to tensile or 
compressive loads of equal magnitude, but opposite sign, at its ends. 

3.3. Equilibrium Equations. It follows from (|3.5|) . (j3.6|) . (|3.15|) . and (|3.14|) that the force-based 
quasicontinuum equilibrium equations 

F^^^{r) + fj = 0, j = -N,...,N + l, (3.18) 

are equivalent to 

^Pf{r) = ^j, j = -N-l,...,N+l. 

We recall that ^n+i = since we have assumed that the external forces satisfy the condition of no 
resultant force (j3.12p . Therefore, if r satisfies the force-based quasicontinuum equilibrium equations 
([3T8]) . then ipN+i{r) = 0, and by ([MD we have Vf (r) = ipf (r) for all j = -A^ - 1, . . . , iV + 1. 
Thus, we see that if r satisfies the force-based equilibrium equations (j3.18p . then 

V'f(r) = $,-, j = -N-l,...,N+l. 

We also recall that if r is symmetric, then 

V^f(r)=V'f(r), j = -N-l,...,N+l. 

Hence, if r is a symmetric solution of 

^f(r) = cl>^-, j = -N-l,...,N + l, 

then r is a solution of the force-based quasicontinuum equilibrium equations (j3.18p . 
If we sum the ghost force iteration equations 

^QCi?(^n+l) ^ pG^^n^ + fj = 0, j = -N,...,N + l, (3.19) 

then we also get the following equivalent corresponding iterative method in terms of the internal 
and external conjugate forces 

V'f (r"+i) + (r") = j = -N-l,...,N + l. 

We also have that if the sequence {r"} satisfies the iteration equations (|3.19p . then the sequence 
{r"} satisfies 

V'f (r"+i) + V'f (r") = j = -N-l,...,N + l. (3.20) 

We finally note that if the sequence {r"} is a symmetric solution of (j3.20p . then the sequence {r"} 
satisfies the ghost force iteration equations (j3.19p . 
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4. Existence of Solutions to the Force-Based Quasicontinuum System 

In this section, we will give conditions on $ = (^^at, . . . , ^n) G such that there exists a 

unique solution r to 

*^(r) = $ (4.1) 

in a domain ^1 C R^^^^, where '^'■^(r) = (^ip^]y{r), . . . ,t/jf^{r)^ G ]R^^+^ has the symmetric form 

given in (|3.10p . We note that we ignore i>^jq_i, V'^+i, *^-Af-i, ^at+i in the above formulation 
since = V'^+i = $-Ar-i = ^at+i = 0. We conclude this section by showing that the 

unique symmetric solution r to ()4.ip for a symmetric $ is a solution of the force-based equilibrium 
equations (jS.lSp . 

For any uniform lattice spacing a = (a, . . . , o) € M^^"*"^, we have a corresponding uniform '^'■^(a). 
Under appropriate assumptions on the atomistic potential, i^(r), and the interatomic spacing, a, 
we will give an explicit neighborhood around a in which we have a unique solution to (j4.ip for any 
external potential, ^, in an explicit neighborhood of ^'^(a). We will show that these assumptions 
on the atomistic potential, 4'{r), are satisfied by the Lennard- Jones potential 

^^r) = jr2-^.- (4-2) 

4.1. Assumptions on the Atomistic Potential, (j){r). We now give the assumptions on the 
atomistic potential, i?i>(r), that are required for our analysis. We will assume that (f){r) € ((0, oo)), 
and we recall that r/(r) = 4''{r) and (j)[r) = (p{r) + 0(2r). We define 

fi{r) = (/>'(r), 

so 

rj{r) = <j)'{r) + 2(/)'(2r) = r?(r) + 2r/(2r). 

We assume that the atomistic potential, 0(r), satisfies the following properties that are graphi- 
cally displayed in Figures [6] and [71 



r/'(r) > for < r < fi and rj'{r) < for r > fi. 


(4.3) 


r/"(r) < for < r < f2 and ri"(r) > for r > f2, 


(4.4) 


fi{r) < for < r < Oo and fj{r) > for r > ao. 


(4.5) 


r/'(r) > for < r < ai and fi'{r) < for r > ai. 


(4.6) 


< ao < fi < f2 < 2ao, 


(4.7) 




(4.8) 



We note that it follows from (|43]) that Fj'^ (ao) = for j = -N, . . . , AT + 1 and *'^(ao) = 

for ao = (ao,...,ao) GM^^+i. 

The local quasicontinuum approximation is simply the case where all of the representative atoms 
are continuum, giving a decoupled system 

*^(r) = 

where ^^{r) = (■0f;^(r), . . . , '^^(r)) is defined by 

V^f (r) = r/(ri), i = -N,...,N. (4.9) 

We thus have that 
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Figure 6. The Lennard- Jones potential (|4.2p demonstrates the prototypical behav- 
ior of (p{r) and its derivatives, rj{r) = (p'{r) and r]'{r) = 4>"{r). 




and we see that the local quasicontinuum approximation is unstable when > ai for some i = 
-N,...,N. 

4.2. Existence and Uniqueness by the Inverse Function Theorem. While the main theorem 
of this section gives explicit conditions on $ for which ^'^(r) = $ is solvable, we begin by showing 
that the inverse function theorem [4] can be used to show that ^^(r) is bijective in a neighborhood 
of r = aQ. Therefore, we must analyze the invertibility of 

' [r?'(r,) + 4r?'(2ri)]<5,,, -N<i< -K, 

[r]'{ri) + r]'{ri + r^.i) + r/'(ri + ri+i)]6ij 
+7]'{ri + rj_i)(5i_ij + 7]'{ri + ri+i)5i+ij 
+ [AT]'{2rK) - iivK + TK-i) - iirx + rK+i)\5Kj 

-7]'{rK + rK-i)6K-ij - ri'{rK + rK+i)5K+ij, -K + 1 < i < K - 1, 

y[n'{n) + Ar^'{2n)]5^,, k<i<n. 



(4.10) 
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Lemma 4.1. //r/'(ao) +8r/'(2ao) > 0, then ^^(r) is bijective in a neighborhood of slq. 

Proof. We will show that D^^(ao) (where {D'^^)ij = Djipf) is nonsingular by demonstrating 
that it is strictly diagonal dominant [20] (with positive diagonal), that is, 

DS!{si(^) > I^J-V-f (ao)I, i = -N,...,N. 

Looking at the rows of D^^(ao), all of the nearest-neighbor terms are on the diagonal. If we 
collect all of the next-neighbor terms on the right hand side of the above inequality, we find that 
it is sufficient to show that r]'{ao) > 4|?7'(2ao)| for rows i = —N, . . . , —K and i = K,. . . ,N and 
r]'{ao) > 8|r/'(2ao)| for rows i = —K + 1, . . . ,K — \. Since 2ao > n by assumption (14. 7p . we have 
that r/'(2ao) < by (j4.3p . Thus, the condition r]'{ao) + 8r]'{2ao) > implies the strict diagonal 
dominance of ()4.10p . Therefore, we have that Z)^^(ao) is nonsingular, and so by the inverse 
function theorem there exists a neighborhood of slq in which ^^(r) is bijective [1]. □ 

4.3. Existence and Uniqueness for General External Forces. In the theorems of this sec- 
tion and the next, we will use the continuation method [1] to find explicit conditions under which 
^■'^(r) = $ has a unique solution. The idea of the continuation method is to start with the uncou- 
pled local quasicontinuum system, construct a homotopy transformation from the local quasicon- 
tinuum system to the force-based quasicontinuum system, and show that existence and uniqueness 
persists through the transformation. We use the following well-known lemma to show that unique 
solvability of the local quasicontinuum system implies existence and uniqueness of a solution to the 
force-based quasicontinuum system [4]. 

Lemma 4.2. Let U C R^iv+i an open, bounded set. Suppose that f, g G C"'^(il; M^^"^^) and the 
homotopy h(r,t) = (1 — t)f(r) + tg{r) for t G [0, 1] satisfies 

(1) f(r) = has a unique solution in $7, 

(2) det Z)rh(r, t) / in J7 x [0, 1], 

(3) h(r, t) / for every (r, t) € x [0, 1]. 

Then there is a unique solution r G satisfying g(r) = 0. 

Proof. Since detL'rh(r, t) ^ for all (r,t) E x [0,1], we know by the implicit function theo- 
rem [4] that the solution of h{r,t) = in a neighborhood of a particular solution (r*,t^,) satisfying 
h{r^,t^) = can be written as a function r(t) in a neighborhood of (r*,t^,). Thus, by compact- 
ness, the curve can be continued until it leaves the region x [0, 1]. Since we have assumed that 
h(r,t) 7^ on di} x [0, 1], there is a one-to-one correspondence between solutions of f(r) = and 
solutions of g(r) = 0. Therefore, the existence of a unique solution of f(r) = for r G implies 
that g(r) = has a unique solution for r G il. □ 

We define the homotopy 

h(r,t) = (l-t)[*^(r)-*]+t[*^(r)-*], tG [0,1], 

where ^^(r) is the local quasicontinuum system ()4.9p . The next lemma uses the properties of the 
local quasicontinuum system to simplify the hypothesis for the previous lemma. 

Lemma 4.3. Let r^ and rjj satisfy < r^ < ru < ai, and set U = (r/,, rj/)^^+-'^. Suppose that for 
i = —N, . . . ,N we have that 

hi{r,t) > if ri = rjj and r G dO., t G [0, 1], 

hi{T,t) <Oifri = rL and r G (9il, t G [0,1]. ^^'^^^ 
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If D^^(r) is strictly diagonally dominant for all r € Jl, then there exists a unique solution to 
*^(r) = 

Proof. We show that the above conditions are enough to satisfy the hypotheses of Lemma 14.21 We 
first show that '^'■^(r) = $ has a unique solution in $7. Since the system is decoupled, we need only 
demonstrate that the scalar equations 

7?(r,)+2r?(2r, ) = j = -N,...,N, (4.12) 

have a unique solution for < rj < rij. We have from the hypothesis ()4.1ip on /ij(r,t) at t = 
that r]{rL) + 2r]{2rL) -^j <0 and r]{ru) + 2r]{2ru) - > 0. Hence, a solution rj to (I4.12|) exists 
by the intermediate value theorem. Since ??(r) = r]{r) + 2r]{2r) is increasing for < r < rjj < ai 
by (|4.6p . the solution to ()4.12p must be unique. 

The hypothesis (14. lip implies that h(r, t) ^ for every (r,t) E x [0,1]. Finally, to show 
that det Dj.h{r, t) ^ for all (r, t) G x [0, 1] it is sufficient to demonstrate the strict diagonal 
dominance of Drh(r, t) for all (r, t) G x [0, 1]. 

Now 

Djhi{v,t) = (1 - t)I),#(r) + tDjijfir), 

and 

D^i;f{r)=f,'{ri)5i,. 

Since < rj < ai, we have by (j4.6p that fj^ri) > 0. Thus, we have that Djil)f{r) is strictly diag- 
onal dominant (with positive diagonal). We can then conclude that Djhi{r,t) is strictly diagonal 
dominant (with positive diagonal) if Djipf (r) is strictly diagonal dominant (with positive diagonal) 
since Djhi{r, t) is then the sum of two strictly diagonal dominant matrices (with positive diagonal). 
Thus, by Lemma 14.21 the equation ^^(r) — $ = has a unique solution in i7. □ 

We next turn to giving results that allow the calculation of an explicit neighborhood of a for 
which '^'■'^(r) is bijective. We first give a condition on < r^ < rjj such that D4/^{r) is strictly 
diagonally dominant for all r G where Q = {r^jru)"^^^^ . 

Lemma 4.4. Suppose that r^ and rjj satisfy 

— <rL < ru, 

ri'{ru) + l2ri'{2rL) > 0. (4.13) 

Then ru < ai, and L'^'^(r) is strictly diagonally dominant for all r £ Q where $7 = {rL,ru)'^^~^^ . 

Proof. First, we note that since r^ > (|4.4p implies that 12rj'{2rL) < 12r]' {2ru). Also, (j4.4p gives 
that the next-neighbor terms are negative, so that rj'{ru) + 4r]'{2rij) > r]'{rjj) + 12r]'[2ru) > 0. 
Therefore, by (14. 7p we can conclude that ru < ai. 

We can obtain from (14.10p by summing all of the next-nearest-neighbor terms in each row that 
L'^'^(r) is strictly diagonally dominant if 

r]'{r) > I2\ri'{2s)\ for ah r, s G {ri^ru). 

(We note that the hypothesis of Lemma l4. 1 1 required the weaker condition ri'{aQ) > 8|r/'(2ao)| since 
we were able to utilize the cancellation of terms in the expression 

Dxi^fir) = Ani'{2rK) - r]\rK + rx-i) - r]'{rK + rx+i), i = -K + 3, . . . - 3, 

when it was evaluated at r = ao to obtain that Dxipf (ap) = 2r/'(2ao) for i = —K + 3, . . . , — 3.) 
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We have by (gS]) and (23]) that 

??'(r) > and r]'{r) is decreasmg for r < fi, (4-14) 

and we have that 

/ ^2 

f] (2s) < and r]'{2s) is increasing for s > —. (4-15) 
It follows from (j4.15p that to prove strict diagonal dominance it is sufficient to show that 

7]'{r) + 12rj'{2s) > for all r, s S {rLifu)- 
We have from (I4.14|) and ()4.15p that the above condition follows from the hypothesis (14.130 . □ 
Theorem 4.1. Suppose that ri and ru satisfy 

y < r-L < ru, (4.16) 
r]'{ru) + 12r?'(2ri) > 0. (4.17) 

// 

7/(rL) + 47/(2rL) - 2??(2rt/) < < ri{ru) + 4r/(2r[/) - 27?(2rL), j = -N,...,N, 
then ^^(r) = $ has a unique solution r in = {rL,ru)'^^~^^ ■ 

Proof. We have by Lemma 14.41 that L'^'^(r) is strictly diagonally dominant for all r £ Q,. By 
Lemma 14.2^ we need only show that 

hi{r,t) > a ri = ru and r G dQ, t € [0,1], 
hi{r,t) < if ri = ru and r e dil., t G [0,1]. 

If we look at the rows of ^^(r), we see that there is always one nearest-neighbor term and at most 
four positive and two negative next-nearest-neighbor terms (j3.10p . We also note that 2rL > f2 > fi 
by (|4.16p and (|4.7p . so we have by (j4.3p that r/(2r) is deceasing for r > r/, and 

27?(2ri) - 27]{2ru) > 0. 

We can thus estimate h(r, t) on the boundary to get 

min h,{r,t)= min - t)[ipf {r) - <^>j] + t[i;f {r) - <^j] 
rGail,rj=ru rGail,rj=ru 

>{l-t) [r]{ru) + 27]{2ru) - <l>j] + t [rj{ru) + 4r/(2rc/) - 2j?(2rL) - ^j] 

= {l-t) [2r]{2rL) - 2r,{2ru)] + [viw) + 4r/(2rc/) - 27]{2rL)] - > 0, , , 

(4-18) 

max /ij(r, t) = max (1 - t)[^'f{r) - <^j] + t[^^ (r) - <i>j] 

r£dQ,rj=r]^ r£d^l,rj=r]^ 

<{l-t) [r,{rL) + 2rj{2rL) - ^j] + t [7]{rL) + 4??(2ri) - 2r,{2ru) - ^j] 
= (1 - t)[27]{2ru) - 2r]{2rL)] + [r/(ri) + 47/(2ri) - 2r]{2ru)] - < 0. 

Therefore, by Lemma 14.31 we have a unique solution to ^^(r) = $ in $7. □ 

We now consider the case in which the external forces, fj, are anti-symmetric (I3.16p . so that the 
external potential $j is symmetric (j3.17p . We then have that the solution r to (r) = $ obtained 
in the proof of Theorem 14.11 is symmetric, and will also solve the equilibrium equations (j3.18p . 
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Corollary 4.1. Suppose that and rjj satisfy 

Y < r-L < ru, 
ri'iru) + l2ri'{2rL) > 0, 
and that the external forces, fj, are anti-symmetric (j3.16p . // 

riiri) + 47?(2rL) - 2r]{2ru) < < v{ru) + 4r/(2rc/) - 2ri{2rL), j = -N,...,N, 
then the equilibrium equations 

F^^^{r) + fj = 0, j = -N,...,N + l, 
have a unique symmetric solution r in U = {rL,ru)'^^~^^ . 

Proof. We first note that since the external forces, /j, are anti-symmetric, the external potential $ 
is symmetric. We consider the solution r{t) to the homotopy continuation 

h (v{t),t) = (1 - {v(t)) - *] + (r(t)) - t G [0, 1], (4.19) 

that we have from Lemma |4.2[ The unique solution to the decoupled local quasicontinuum system 
^^(r) = $ is symmetric, therefore r(0) is symmetric. We obtain by differentiating (|4.19|) that 

Vrh (r(t), t) n (t) + ht {r{t),t) =0, t G [0, 1], 

where 

Vrh (r(t), t) = (1 - t)Vr*^ (r(t)) + tVr*^ (r(t)) 

and 

ht (r(t),t) = *^(r(t))-*^(r(t)). 

All terms above are symmetric whenever r{t) is, and therefore the symmetry of r(0) implies that 
rt{t) and r(t) are symmetric in the whole interval [0, 1], and in particular r(l) is symmetric. The 
proof that r is symmetric is completed by observing that r(l) is the unique solution to ^^(r) = 

Since the unique solution r to ^^(r) = $ is symmetric, we can conclude that r is the unique 
symmetric solution to the force-based equilibrium equations 

F^^'^{r) + fj = 0, j = -N,...,N + l. □ 

We now apply Corollary 14. II to the Lennard- Jones potential (14. 2|) . 

Corollary 4.2. We assume that the external forces, fj, are anti- symmetric (j3.16p . For any ^ < 

ru < ai, let 

^. = maxfef^)*V (4.20) 



16ri'{ru) 



If ri < ru, then the equilibrium equations 

pQCF 



^7"^ (r)+/i = 0, j = -N,...,N+l, 
have a unique symmetric solution r in = {rL,rij)'^^'^^ whenever 

Viri) + 47?(2rL) - 2r]{2ru) < < v{ru) + 4r/(2rc/) - 27?(2rL), j = -N,...,N. 
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Proof. To prove this corollary using Theorem 14. H we need to show that (j4.17p holds. Evaluating 
(j4.17p for the Lennard- Jones potential (j4.2p . we get 



rj'iru) + l2i]'{2rL) = v'{ru) + 12 

63 



156 84 



214^14 28r^j 



Therefore, holds if vl > [jq^) • □ 

After solving for in (|4.20p . we have to check the additional hypothesis (j4.16p that rj^ < rjj, 
as it is not true for every rjj € (r2/2,ai). However, this is not a very restrictive assumption since 
it is true for all rjj < 1.1003, whereas ai = 1.1059. The lower end for the interval is ^ = 0.6085. 

Using Corollarv l4.2l we now find a symmetric region about $ = 0. Solving numerically, if we take 
ri = 0.9700 and ru = 1.0883, we can conclude that for any $ satisfying —2.62 < $ < 2.62, we can 
uniquely solve the force-based quasicontinuum equilibrium equations (j3.18p . For the Lennard- Jones 
potential 57(01) = 2.781, so this symmetric region extends quite close to the load limit. 

Remark 4.1. The techniques of this section can be applied to the analysis of the fully atomistic 
model or the constrained atomistic quasicontinuum approximation, and the analysis will be sim- 
plified as there are no non-local conjugate forces. In both cases, the continuation from the local 
quasicontinuum approximation can be used. 

We also note that the interfacial terms satisfy the convexity condition 

5,-(r,_i, r,) = -l</.(2r,_i) + (/.(r,_i + r,) - ^<P{2rj) > 0, j = -N + 1, . . . , N, 

for r in the region defined by (j4.16p in Theorem 14.11 and Corollary 14.11 since ()4.15p holds in this 
region. 



5. Convergence of the Ghost Force Iteration 
We now give a similar analysis for the iterative equations ()3.20p which in vector form are 

*^(r"+^) + *^(r") = (5.1) 

where ^^(r) = (^'4^^j^{r), . . . , ip^{r)^ G M^^+^ has the symmetric form given in ()3.4p . and ^''^(r) = 
(^ip'^j^(r), . . . ,^^(r)^ G M^^+^ has the symmetric form given in (13.1ip . We also note that we can 
ignore V'^at-ij V'^^+ii i^'^N-i^ ^^+1 ^'^^ ^-N-i^ ^N+i in the above formulation since 

We will determine Q C and D C R^^+^ such that there is a unique r"+^ G Q, that satis- 

fies (|5.ip whenever r" G and $ G D, and we will show that the induced mapping from r" to r""^^ 
is a contraction mapping. 

First, we need to compute D'^^(r) and D'^'-'{r). The following expressions are rather complex, 
but all we will need to recognize from them is that the number of first and second neighbor terms 
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in each row. We have 

qr?'(r,)+4r/'(2ri)]5,j, 
[r/'(r,) + 4r/(2ri) + ir?'(ri + ri^i)]Sij 

+ [^??'(ri + rj_i)]5j_ij, 
[r]'{ri) + 2r]'{2ri) + ir/'(ri + r,_i) + y{r, + ri+i)]Si 

+ [^ri'{ri + rj„i)]5i_ij + [^r?'(ri + ri+i)]5j+ij, 
[r/'(ri) + ??'(ri + ri_i) + i7/(ri + ri+i)](5ij 

+ W{ri + ri_i)]Ji_ij + [^r/'l^i + ri+i)]Si+ij, 
Wiri) + ??'(rj + ri_i) + r/'(ri + ri+i)]6ij 
+r]'{ri + rj_i)(5i_ij + r]'{ri + rj+i)5j+ij, 



and 



-\'r]'{rK + rK+i)5Kj - ^v'i^K + rK+i)SK+ij, 
[2r]'{2rK) - \'r]'{rK + rx-i) - \i{rK + rR+iil^Kj 

-\r]'{rK + rK-i)5K~ij - ^H^k + rK+i)6K+ij, i = K, 
[V(2r/^) - lri'{rK + rx-i) - r]'{rK + rK+i)]SKj 

-Iv'irK + rK^i)SK-ij - v'ifK + rK+i)5K+ij, i = K - I, 
[Ar]'{2rK) - i{rK + tr-i) - v'i^K + rK+i)]6Kj 

-r]'{rK + rK^i)6K~ij - i{rK + rK+i)6K+ij, 



K + 2 <i < N, 
i = K + l, 
i = K, 
i = K -1, 

-K + 2 <i < K - 2, 
-N <i < -K + 3, 

K + 2 <i < N, 
i = K + l, 



(5.2) 



(5.3) 



K + 2<i<K-2, 
N <i < -K + 3. 



We recall that the maximum norm for r G M^^+^ is given by 



r||oo = , max \ri 



Theorem 5.1. Suppose that and ru satisfy 

— <rL < ru, 



(5.4) 



r]'{ru) + l37]'{2rL) > 0, 
and that $ satisfies 

riiri) + 4r/(2rL) - 2r]{2ru) < < r]{ru) + 4r/(2rc/) - 2r/(2rL), j = -N,...,N. (5.5) 
Then for every r" S = {ri,ruY^^^ there is a unique r"+^ € $7 such that 

*^(r"+i) + *^(r") = *. (5.6) 
We also have that the induced mapping r" r^+i ig ^ contraction and satisfies the inequality 



|„n+l _ „n+l|| ^ 



8\r]'{2rL)\ 



where we have from (|5.4|) that 



r]'{ru) - 5|r/'(2ri)| 
8|V(2ri)| 



(5.7) 



7?'(r[/)-5|r/'(2rL)| 



< 1. 
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Proof. The first part of tfie proof will be very similar to the proofs of Lemma 14.41 and Theorem 14. H 
as we will again satisfy the hypotheses of Lemma [4.31 First, we note that (j5.4p implies that rfj < ai 
by Lemma [4.41 We prove strict diagonal dominance for Dj^f'(r) given by (15. 2p . For this argument, 
we only need to show that r]'{r) > 5\ri'{2s)\ whenever r,s £ {fL^fu)^ or since r]'{2s) < 0, we need 
only show 

r]' (r) + 5r]' {2s) > for r, s G (r^, r^/). 
We will need the factor 13 in the hypothesis (|5.4|) to prove that the mapping is a contraction, 
which is why the hypothesis is as strong as it is. We further note that from (I4.14p . ()4.15p . and the 
hypothesis ()5.5p that we have 

7]'{r) + 5??'(2s) > r?'(r) + 13r/'(2s) 

> v'{ru) + 13r/'(2ri) > 0, r,s € {rL,ru)- 

Thus, we have established that detD^^(r) > in 17. 

We now verify that h(r, t) satisfies the condition (14. lip , that is, it does not vanish on 9$7 x [0, 1]. 
In the proof of Theorem 14.11 we analyzed h(r, t) on dfl x [0, 1] by grouping the second-neighbor 
terms. From this point of view (|4.18p . the iteration problem ^^(^r"'*'^) + ^'^(r") = is identical 
to the problem, ^{r) = Thus, we have that 

min /i,-(r"+i)> min (1 - t)[#(r"+i) - ^>,-] + t[^f (r"+^) + ^^(r") - 

rj=ru rj=ru 



> (1 - t)[2r]{2rL) - 2ri{2ru)] + t [rj{ru) + 4r/(2rc/) - 2rj{2rL) - <^j] > 0, 
max /i,-(r"+^) < max (1 - t)[#(r"+i) - + t[V;f (r"+i) + ^^^(r") - 



< (1 - t)[2ri{2ru) - 2??(2ri)] + t [T]{rL) + 4r/(2rL) - 2r]{2ru) - ^j] < 0. 

Therefore, we have proven that h{r,t) satisfies the condition ()4.1ip . that is, it does not vanish on 
di} X [0, 1]. We have now verified the hypotheses of Lemma [4.3l to conclude the existence of a unique 
solution r"+^ € to the iteration equation (j5.6p . 

To prove that the mapping (15. 6p is a contraction, we suppose r", r"+^, s", s""''^ G Q, satisfy 

*^(s"+i) + *^(s") = 

We then have that 

*^(r"+i) - *^(s"+i) = *^(s") - *^(r"). (5.8) 
By the fundamental theorem of calculus, we have that 

#G(^n) _ ^G(gn) ^ L^(r",s")(r" -s"), 

where 



(5.9) 







L^(r", s") = d4/^ (r" + 0(s" - r")) dO. 
Jo 



We then have by <^M) and <^M) that 

L^(r"+\s"+^)(r"+i -s"+i) = -L^(r", s")(r" - s*^ 
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To show that L^(r"+^, s""*"^) is nonsingular and estimate its inverse, we define the diagonal 
matrix, L^(r"+^, s"'"'"^), that contains the dominant nearest-neighbor terms of L^(r"'+^, s""*"^) by 
^2^D(j.n+i^gn+i)^^^ = Lg(r"+i,s'^+i) where 

-1 



Jo 



n+l 



)) de 



6ij. 



Since r" + 0(s" — r") £ Q, = (r^,, rf/)^^"*"^ if r", s" G and 6 G (0, 1), we can conclude from 
(fil^ and Km\ that 



|L^(r"+i,s"+i)-i| 



< 



r/'(rt/) 
<5|r/'(2ri)|. 



(5.10) 



where the matrix norm that is induced by the ||r||oo vector norm is 

JV 

iLlIno = max ^ \Lij\. 



loo ^ max 

i=-N,...,N 



j=-N 



We have that 



so it follows from (jS.lOp that L'^(r"^^, s"'^^) is nonsingular and we have the estimate 



Hence, can state that 



< 



|L^(r"+i,s"+i; 



-ii 



1 - ||L^(r"+i,s"+i)-i||ooI|i^^(r"+i,s"+i) -L^(r"+i,s"+i 



(5.11) 



(5.12) 



From ()5.3p . we can obtain from ()4.15p the estimate 

||^^(r",s")|U <8|7?'(2rL)|, 
so we have from (j5.12p and (jS.lip that 



„n+l 1 



Jllooll^ * ||oo 



ir- ^ - s-^Hoo < ||L^(r"+\s"-^^)-^L^(r",s-)||^||r- - s 



< ||L^(r"+\s"+^)-^||oo||L^(r",s")||oo||r"-s"||oo 



< 



We have from (15.411 that 



E/^n+l gn+l^-l| 

^IV(2rL)| 

r 

rj'{ru) - 5|?7'(2rL)| 

8W{2rL)\ 
rj'iru) - 5\r]'{2rL)\ 



< 1, 



so the mapping r" — > r"'^'^ is a contraction and 



n+l „n+l 1 



< 



8|r/'(2ri)| 



r]'{ru) - 5\r]'{2rL)\ 



□ 



We can prove the following corollary of Theorem 1 5. II for the ghost force iteration by an argument 
similar to that used to derive Corollarv 14.11 from Theorem 14.11 
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Corollary 5.1. Suppose that the external forces, fj, are anti- symmetric (j3.16p . that r^ and rjj 
satisfy 

— <rL < ru, 

ri'iru) + 137?'(2rL) > 0, (5.13) 

and that ^ satisfies 

Viri) + 47?(2rL) - 2r]{2ru) < < v{ru) + 4r/(2rc/) - 27?(2rL), j = -N,...,N. 

Then for every symmetric r" G = (r/,, rf/)^^"*"^ there is a unique symmetric r"+^ G such 
that 

pQCE ^^n+i^ + if (r'^) + /i = 0, j = -N,...,N + l. 
We also have that the induced mapping r" — > r""*"^ is a contraction and satisfies the inequality 

-s"+^|U < y^':^?' J |r"-s"|U, (5.14) 
r]'{ru) - 5\r]'{2rL)\ 

where we have from (|5.13p that 

8jV(2ri)| 



< 1. 



r]'{ru) - 5|r/'(2rL)| 

The mapping r" — s- r"+^ converges to the unique symmetric r in Q that satisfies the force-based 
quasicontinuum equations 

F^^^{r) + fj = 0, j = -N,...,N + l. 

We now apply Corollary 15.11 to the Lennard-Jones potential as in the previous section. This 
time, we not only need to satisfy the basic inequality (|5.4p on ri and ru, but we also need to verify 
that the contraction constant in (j5.14p is less than 1. So, we pick 7 G (0, 1) and solve for 

8W{2rL)\ 
ri'iru) -5\ri'{2rL)\ ^' 
Using the same argument as in Corollarv 14.21 for any ru G (f2/2, ai), we choose 



ri = max 



'f2 / 156(5 + 8/7 ) 
2 ' V 256r/'(rc/) 



If the resulting r^ is less than r^/, we then have a region for symmetric $ for which the iteration 
is well-defined and a contraction. 

Using the above with contraction constant 7 = ^, we find that for any symmetric $ G (—2.56, 2.56)^^"*"^ 
and symmetric G (.9706, 1.0771)^^+^ there is a unique symmetric r"+^ G (.9706, 1.0771)^^+^ 
that satisfies the ghost force iteration equations 

pQCE^^n+l^ + Ff{v^) + /, = 0, j = -iV, . . . , AT + 1. 

We can finally obtain by taking s" = r, where r is the unique symmetric solution to the force-based 
quasicontinuum equations, 

i^?''^(r) + = 0, j = -iV,...,iV + l, 



that 

llr"+^ -rlU < -llr" 
2 
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